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Abstract 

We derive an analytic expression for the shear viscosity of an ultra-relativistic gas in presence 
of both elastic 2 — )• 2 and inelastic 2 3 processes with isotropic differential cross sections. 
The derivation is based on the entropy principle and Grad's approximation for the off-equilibrium 
distribution function. The obtained formula relates the shear viscosity coefficient 77 to the total cross 
. sections o"22 and CJ23 of the elastic resp. inelastic processes. The values of shear viscosity extracted 

using the Green-Kubo formula from kinetic transport calculations are shown to be in excellent 
agreement with the analytic results which demonstrates the validity of the derived formula. 
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I. INTRODUCTION 



Dynamics of the hot and dense state of nuclear matter, produced in heavy-ion colhsions 
at modern accelerator facilities like RHIC at Brookhaven National Laboratory or LHC at 
CERN, can be described by means of kinetic transport theory or relativistic hydrodynamics. 
In fact, these two approaches have been widely used for research of the properties of the 
so-called Quark-Gluon Plasma (QGP), early-time dynamics of which determines the prop- 
erties of hadronic and electromagnetic secondaries measured experimentally [liJal. Whereas 
the kinetic transport theory is based on a microscopic paradigm and describes evolution 
of the phase-space distribution function by solving the corresponding transport equations, 
hydrodynamics is a macroscopic theory describing space and time evolution of macroscopic 
fields like the energy-momentum tensor and particle number four-vector. For the hydro- 
dynamic theory the transport coefficients are an extrinsic input and must be calculated 
from the underlying field theory, which for QGP is the Quantum Chromodynamics (QCD). 
Calculations of the transport coefficients are based on the formal correspondence between 
the kinetic theory and the hydrodynamics, which can be obtained by integration of the 
microscopic equations over the phase space. For example, such correspondence is discussed 
in {4-8]. Moreover, kinetic transport calculations have often been used as a benchmark 
for hydrodynamic models, as for instance in Refs. 9l-ll l| . which again requires a consistent 
connection between the cross section on the one hand and transport coefficients on the other. 

For a strongly interacting QGP, the shear viscosity coefficient has been calculated in the 



perturbative regime as function of the coupling parameter as for instance in Refs. 12 1 



13 



15|. Such calculations account for the radiative processes occurring in the QGP, which have 
an important contribution to the isotropization of the medium. The role of inelastic processes 
was understood by means of the transport cross section or transport rate [l^, which are 
analogous to the total cross section and collision rate but are additionally weighted by the 



angle of outgoing particles. In contrast, calculations reported in Refs. [16| do not see a 
strong contributions of the inelastic 2 — )■ 3 processes. In Ref. [l^ the shear viscosity of 
a QGP was calculated employing two-particle scatterings, i.e. without taking into account 
the radiative processes. 

In this work we derive an expression which relates the shear viscosity coefficient of an 
ultrarelativistic gas to the total elastic and inelastic cross sections in the medium. For our 



derivation we consider isotropic, i.e. momentum- independent, matrix elements for the coUi- 
sional processes, thus bypassing the conceptual difficulties connected with the momentum- 
dependent matrix elements in a (perturbative) QCD medium. To our best knowledge, ana- 
lytic derivation of the shear viscosity for 2 ^ 3 inelastic processes with isotropic differential 
cross sections is completely new. 

This paper has the following structure. In Section [ITl we start with the derivation of 
an expression for the shear viscosity coefficient, for which we use the entropy production 
in a system off equilibrium. In order to obtain a closed analytic expression for the shear 
viscosity coefficient in terms of the involved total cross sections, kinematic integrals will be 
evaluated. This is discussed in Section IIIH where the final expression is also presented. In 
Section |IV] the results of the analytic formula are compared with Green-Kubo based kinetic 
transport calculations of the shear viscosity coefficient. Finally, Section |V] contains summary 
and conclusions. 



II. DERIVATION OF THE EXPRESSION FOR THE SHEAR VISCOSITY COEF- 
FICIENT 



The derivation we follow is built upon the entropy principle, i.e. the second law of 
thermodynamics is used. In presence of dissipative fields, of which we will consider only 
the shear stress tensor throughout this work, the entropy production (?^s^ must have a 



non- negative form, and hence can be written as discussed e.g. in Ref. 19| : 



a .A« = !l>!l_ + jlnA>0. (1) 

In the latter expression ir^'^ denotes the shear tensor and rj the shear viscosity coefficient; 
T is the temperature and A the fugacity, which are both related to the chemical potential, 
InA = ^. J denotes a source for particle production and annihilation and vanishes in the 
situation where the net particle number does not change. 

Note that rj in Eq. ([1]) is identified with the shear viscosity since one can show that in the 
lowest-order, i.e. Navier-Stokes, dissipative hydrodynamic theory, where n^" = 2r]V^^u^'^, 
the entropy production indeed takes the form ([1]). Since higher-order theories (such as Israel- 
Stewart [J]) are supposed to converge to the Navier-Stokes form, the entropy production is 
always written in the form ([T]) and a consistent dynamic equation for tc^'^ can be found, like 



flU 



e.g. in Ref. 

On the other hand, the entropy production can be obtained from the total divergence of 
the entropy current defined in the kinetic theory as 



p^f{x,p) [ln/(a;,p)-l] dT 



(2) 



where f{x,p) denotes the phase-space distribution function and dT = d^p/{2TT)^/E. For 
systems away from kinetic equihbrium the distribution can be written introducing a small 
deviation (f){x,p) from the equilibrium form fQ{x,p): 



f{x,p) = fo{x,p) ■ (1 + (p{x,p)) 



(3) 



Some constrains must be applied to Eq. in order to preserve consistency. For instance, 
if moments of f{x,p) from (j3]) are calculated, the particle number four-vector A^'^ and the 
energy-momentum tensor T^'^ must be recovered. This means that the deviation (j){x,p) 
must depend on the deviations ST'^'^ and 5A^^ of the energy-momentum tensor resp. particle 
number four-vector from their ideal forms. The explicit functional dependence is given by 



the following expression, as obtained in scope of Grad's method of moments 



= c(e,p,T) ■ Ti^uP^p'' ^ 
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22|: 



(4) 



Note that the momentum-dependence of (j) can be extended beyond the second power, as 



demonstrated in Ref. 



24j |. However, our derivation is restricted to the original approach of 



Grad, for which the series must be truncated at second power of momentum. In the kinetic 
theory, dynamics of f{x,p) is governed by a transport equation, such as the Boltzmann 
Equation 

p''d,f{x,p)=C[f{x,p)] , (5) 
with the functional C [f{x,p)] on the right hand side denoting the collision term. Combin- 



ing Eqs. (E 
coefficient [13|, 



22| 



and ([T]) one obtains the following expression for the shear viscosity 



III/ 



2c(e,p,T) ■ Ttt^^P 



with V^" denoting the second moment of the collision term 



flU 



pVC[/]rfr 



(6) 



(7) 



An obvious property of the obtained expression is the exphcit dependence on the shear tensor 
components. But since the shear viscosity coefficient is a medium property and should not 
depend on the dynamics encoded in tt'^", Eq. ([6]) should be considered in the limit tt'^'^ — )■ 0. 

The functional V^" contains information about the coUisional processes in the medium. 
We consider only local interactions between constituents of the medium and assume them 
to be Boltzmann-particles, i.e. quantum statistics effects are neglected. If both elastic 
and inelastic interactions are considered, the second moment of the collision term takes the 
following form: 



^ j dwidw2^ j dw[dw'2PiPif[f'2\Mv2'^i2?{2TiY5^^\p\+p'2-pi-p2) 



^ j dwidw2^ j dw[dw'2PiPlfif2\Mi2^i>2'\'^{27ryS^^\p[+p2-pi-p2) + 
^ j dwidw2dw3^ j c?wirfw2PiPi/{/2l^i'2'^i23p(27r)^5(^)(p'i +P2 -Pi -P2 -Ps) + 
+ j dwidw2^ j dw[dw'2dw'^p^p\ f[f2f[i\Mv2'w ^12? {'2t^Y5^^\p\+p'2+p'z- Pi - P2) - 

~\ j ^^1^^2rfW3^ j dw[dW2p'lpl f I f2f 3] Mi23^1' 2' C^Try S^'^^pi + P2 + P3 - P[ ' p'2) " 

dwidw2^ j dw[dw'2dw'^p'lp\fif2\Mi2^i,2'3'?{'2T^Y5^^\pi+P2-p'i -p'2 -P'3) 



(8) 

with dw = |c/r = 2(inf^- '^^^ ^^^^ terms in ([8]) account for elastic 2 — )• 2 processes, 
whereas the last four terms capture the inelastic 2 — 3 and reverse 3 — )■ 2 processes. The 
numerical factors 1/2 and 1/6 account for multiple counting of identical Boltzmann particles. 

In order to evaluate the integrals in ([H]) analytically, the distribution functions will be 
replaced by the near-equilibrium approximations and only terms up to ffist order in 
will be considered: 

" " " (9) 



i \ j / i 



The equilibrium distribution /q is the Boltzmann distribution 

/o = d.e-P'^"°/^+^/^, (10) 

with the degeneracy factor d. 
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Analytic evaluation of the integrals in (jS]) can be simplified if a specific symmetry for 
the geometry of the system is considered. Since the shear viscosity coefficient is a material 
property, its value cannot depend on the form of the energy-momentum tensor distortion, 
i.e. any plausible form for the shear tensor can be considered to evaluate (j8]). The simplest 
geometry that can be considered for the shear tensor, which has to be traceless, is a diagonal 
from with one independent component vr: 

vr^,, = cim5f(0,7r/2,7r/2, -tt) . (11) 

This form of the shear tensor implies deviation from isotropy in z-direction and isotropy in 
';he transverse plane. Systems with this kind of symmetry have been studied e.g. in Refs. 
Ill , lol . With the shear tensor given by Eq. fllip . the expression for the shear viscosity 
(E]) as well as the collision term moment (|H]) are significantly simplified, since now vr is the 
only unknown parameter, for which the vr — )■ limit must be considered. In fact, with the 
introduced form of the shear pressure tensor the correction (j4]) to the distribution function 
becomes (for sake of compactness we will suppress the dependence of c(-) on e,p,T in the 
notation) 

(f){x, p) = c ■ TT ■ - . (12) 

Now the contraction tt^^V^'^ in can be rewritten as follows 



^ f ■ ' 'I 1.1 / 1 2 2 \ H H 



- / dwidw2dw^dw2 \i^Pi^T-Pi,z ] ■ /o,i/o,2 X 



X (l + C • TT . (ip'iV - P'L) + C ■ ^ • i\p'lT - Pt)) ■ |M1'2'^12P(2^)^5(^) (..) 

^ dwidw2dw[dw2 \ \pi.T - Pi A ' /o, 1/0,2 x 



X (1 + c • TT . {\pI^ - pL) + c • ^ • {\pIt - P2, j) ■ |Mi2^i'2f (2^)^<5(^) (..) + 



^ J dwidw2dw3dw[dw'2 {^pl,T-Pi.zj /o,i/a2 x 
X (1 + c • vr . - P'l^) + c • ^ • ilp'i,T - pV) ■ I Afi'2'^i23p(27rf^W (..) + 

^ j dwidw2dw[dw'2dw'3 {]^p1^t - Pi,^ • /o, 1/0,2/0,3 x 

X (1 + c • TT . [\p'l^ - p'l) + c . ^ . {\p'l^ - p'l:) + c • ^ • i}-p'l^ p'l,)\^ ■ |Mi,2'3'^i2p(2^)^<5(^)(. 
^ J dwidw2dw3dw[dw'2 (^^pj^r - Pi,zj • /o, 1/0, 2/0, 3 x 



^ y dwidw2dw[dw'2dw'^ Qpi,T "Pm^ • /o, 1/0,2 x 

(13) 

The parts of the integrals that do not contain c (i.e. those parts where the integration goes 
over /o and not /o ■ 0) have to vanish, since there is no entropy production in an equihbrated 
medium. The remaining integrals will then be of order vr^, just as the expression in the 
numerator of ([6]), so that the shear viscosity coefficient will not depend on vr, as it should 
be. 

III. EVALUATION OF THE KINEMATIC INTEGRALS 

In general the matrix elements |M„_>m| will have a highly non-trivial dependence on 
momenta of the colliding particles and analytic evaluation might be difficult. Here we 
will consider angle-independent matrix elements to obtain the shear viscosity as function of 
total cross section. The total cross sections (J22 and o"23 for elastic resp. inelastic two-particle 



processes can be defined according to l23l |: 

CX22 = J dwi y rfw2|Mi/2'^i2p(27r)^5(pi +P2 - pi -P2) , (14) 



0-23 = ^^ J dwi j dw2 J dw 3] Mr 2' ^123]'^ i2TcyS{pi+p2 + P3-p'i-P2) ■ (15) 

For these definitions the center-of-mass energy squared of the colliding particles s = (^1+^2)^ 
has been introduced. For a process with 3 particles in the initial state a classical cross section 

n 

cannot be defined, but one can define an analogous quantity [2[ 

2^32 = ^ / J '^^2|Mi/2'3'->i2r(27r)''5 {pi +P2 - p[ - P2 - P3) . (16) 

The squared matrix elements |M2^.3p and |M3_j,2p depend on each other by virtue of the 
detailed balance requirement: 

|M2^3|' = rf-|M3^2|'. (17) 

If the matrix elements are angle- independent, i.e. if the scattering is isotropic, as consid- 
ered for the following derivations, one can obtain simple relations between them, the cross 



sections (T22, (^23 and the quantity X32 [2^ by performing the momentum integration in the 
definitions f[T^ - flTB]) . Note that the right-hand-sides of Eqs. f ll4l) - f[T^ are Lorentz- 
invariant and thus can be conveniently evaluated in the center-of-mass frame of colliding 
particle pair: 

|M2^2|' = 327r-s-a22, (18) 

|M2^3|' = 1927r3.a23, (19) 
192 

2:32 = -^^'^23 . (20) 

Here, again, d is the degeneracy factor of the constituent particles. 

With the definitions f ll4p - fll6p half of the terms (the second, fifth and sixth terms) in ([S]) 
are significantly simplified. Calculation of the remaining terms is a bit more involved. The 
first, third and fourth terms in (|8]) contain integrals over the n = 2 or 3 final state momenta 
dwi . . . dwn which can be formally written as 

/ dWi... / dWnGipi . . .Pn)\Mm^n\'^S(^ 

Pinitial J^Pfinal). (21) 

m n 

In the latter expression Q denotes a (polynomial) function of the final state momenta. For 
the matrix elements the relations f|T8|) - fl20l) can be used. However it is important to note 
that the expression fl21l) is not anymore Lorentz-invariant due to presence of the function Q. 
At least for the terms describing a two-particle initial state it is still convenient to transform 
the integral into the center-of-mass frame for evaluation and to apply a back-transform 
afterwards. This means rewriting Eq. (12T1) in the following form 

I dwl... J dw^Qip*^^^, . ■.p*TR,n)\M„^^nmJ2Ptniti.l " E^final) • (22) 



where * denotes quantities in the center-of-mass frame, in which EI + E2 = \/s andpl+p^ = 
0. To rewrite fl^Tl) into the form fl2^ Lorentz invariance of dwi and of the 4-dimensional 
Delta-function can be used. Notation Q{pxri ■ ■ -PrRn) means that function Q{- ) acts on 
momenta that are Lorentz-transformed from the center-of-mass frame into the lab frame: 

E;^ = 7(E*+r/3cM) (23) 
rxR = r + ^^?^/3cA/ ■ r + Icm^cmE* . (24) 

1 + ICM 
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with the Lorentz boosts 



ficM 
1cm 



P1+P2 
El + E2 ^ 



(25) 
(26) 



CM 



For the term describing a three-particle initial state (the fourth term in ([8]) resp. (|T3i) ) 
such transformation into and then back from the center-of-mass frame does not induce any 
significant simplification of the calculation so that the integration must be evaluated directly 
in the lab frame. 

After evaluating all the integrals we finally obtain the following expression for the shear 
viscosity coefficient: 



6 



T 



^ = 7 



5 ^22 + §0-23 + ^ACT23 



(27) 



The second and the third terms in the denominator of the latter equation account for 
contributions of 2 — )■ 3 resp. 3 — ?■ 2 processes. The factor A in the 3 — 2 term denotes the 
fugacity, which is a measure for the degree of saturation of the particle density in a system 



A = n/n, 



eq ) 



(28) 



with = d ■ J e "f^^Z-^ci^p. Setting A = 1 (or, accordingly, yU = in Eq. ( ITU]) ), i.e. 



considering a chemically equilibrated system, we obtain 



6 



T 



^ = 7 



5 cr22 + ^cr23 



(29) 



Neglecting the contribution from inelastic scatterings, this formula is identical with the 
already known results from Refs. 21j and 24i\ . 



IV. COMPARISON WITH KINETIC TRANSPORT CALCULATIONS 

In order to confirm the validity of Eq. fl2U]) . one can compare it to the values of shear 
viscosity extracted from the Boltzmann Approach to Multi-Parton Scatterings (BAMPS) 
using the Green-Kubo formalism. In Ref. we have recently reported, how 

the Green-Kubo formalism can be applied to extract shear viscosity from BAMPS kinetic 
transport simulations. In Fig. [T] we demonstrate the shear viscosity values extracted this 
way from BAMPS calculations for which only inelastic 2 — )■ 3 and 3 — )■ 2 processes were 

9 
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FIG. 1. Shear viscosity in presence of inelastic processes only, extracted from BAMPS using Green- 
Kubo formalism (symbols) and calculated using Eq. (|29|) . The results are obtained employing 
angle-independent (isotropic) differential cross sections. 



considered. The shear viscosity is shown as function of the inelastic cross section (T23. We 
observe a very good agreement of the results extracted from BAMPS (symbols) with the 
results obtained using Eq. fl29|l (lines) derived in this work. We also would like to mention 
that Eq. (129|) confirms the values of the shear viscosity which we obtained in Ref. [3] 
by numerical evaluation of the Equation ([6]). Equation fl29|) implies that if the total cross 
sections are equal, relaxation of the medium towards equilibrium will proceed 1.5 times faster 
if only inelastic processes are involved as compared to the case when only elastic processes 
are involved (since the hydrodynamic relaxation time is directly proportional to the shear 
viscosity coefficient - comp. 5|, llO|). This reflects the fact that the inelastic processes are 
more efficient in driving the medium towards equilibrium, as was demonstrated in Refs. 



13 



18 



19|. 



V. CONCLUSIONS 



We have derived an analytic expression relating the shear viscosity coefficient of an ultra- 
relativistic gas to the total elastic and inelastic cross sections of isotropic scattering processes 

10 



in the medium. From the obtained expression one recognizes that, given all total cross 
sections are equal, isotropic inelastic processes are 1.5 times more efficient than the elastic 
ones in restoring kinetic equilibrium in a system. The values of the shear viscosity obtained 
by Eq. (129|) are in very good agreement with the results of Green-Kubo based calculations 
performed using the kinetic transport model BAMPS. The correspondence between shear 
viscosity and total cross sections of elastic and inelastic processes reported is another step to 
establish a connection between the kinetic transport and dissipative hydrodynamic models. 

ACKNOWLEDGEMENTS 

This work was supported by the Helmholtz International Center for FAIR within the 
framework of the LOEWE program launched by the State of Hesse. BAMPS calculations 
were performed using the resources of the Center for Scientific Computing of the Goethe 
University Frankfurt. 



[1] H. Song and U. W. Heinz, Phys. Rev. C77, 064901 (2008) [arXiv:0712.3715 [nucl-th]]; 

M. Luzum and P. Romatschke, Phys. Rev. C 78, 034915 (2008) [Erratum-ibid. C 79, 039903 

(2009)] [arXiv:0804.4015 [nucl-th]]; 

D. A. Teaney, [arXiv:0905.2433 [nucl-th]]; 

K. Dusling and S. Lin, Nucl. Phys. A 809, 246 (2008) [arXiv:0803.1262 [nucl-th]]. 
[2] Z. Xu, C. Greiner, Phys. Rev. C71, 064901 (2005). [arXiv:0406278 [hep-ph]]; 
[3] Z. Xu, C. Greiner and H. Stocker, Phys. Rev. Lett. 101, 082302 (2008) [arXiv:0711.0961 

[nucl-th]]. 

[4] W. Israel, Ann. Phys. (N.Y.) 100 (1976) 310; J.M. Stewart, Proc. R. Soc. London , Ser. A 357 

(1977) 59; W. Israel, M. Stewart, Ann. Phys. (N.Y.) 118 (1979) 341. 
[5] A. Muronga, Phys. Rev. C69, 034903 (2004). [arXiv: 0309055 [nucl-th]]. 
[6] R. Baier and P. Romatschke, Eur. Phys. J. C 51, 677 (2007) [nucl-th/0610108]. 
[7] B. Betz, D. Henkel and D. H. Rischke, J. Phys. G G 36, 064029 (2009). 
[8] G. S. Denicol, T. Koide, D. H. Rischke, Phys. Rev. Lett. 105, 162501 (2010). [arXiv: 1004.5013 

[nucl-th]]. 

11 



[9] P. Huovinen and D. Molnar, Phys. Rev. C 79, 014906 (2009) [arXiv:0808.0953 [nucl-th]]. 
[10] A. El, Z. Xu and C. Greiner, Phys. Rev. C 81, 041901 (2010) [arXiv:0907.4500 [hep-ph]]. 
[11] I. Bouras, E. Molnar, H. Niemi, Z. Xu, A. El, O. Fochler, C. Greiner and D. H. Rischke, Phys. 

Rev. C 82, 024910 (2010) [arXiv: 1006.0387 [hcp-ph]]. 
[12] P. B. Arnold, G. D. Moore and L. G. Yaffe, JHEP 0011, 001 (2000) [hep-ph/0010177]. 
[13] A. El, A. Muronga, Z. Xu and C. Greiner, Phys. Rev. C 79, 044914 (2009) [arXiv:0812.2762 

[hep-ph]]. 

[14] C. Wcsp, A. El, F. Reining, Z. Xu, I. Bouras and C. Greiner, Phys. Rev. C 84, 054911 (2011) 

[arXiv: 1106.4306 [hcp-ph]]. 
[15] F. Reining, I. Bouras, A. El, C. Wesp, Z. Xu and C. Greiner, Phys. Rev. E 85, 026302 (2012) 

[arXiv:1106.4210 [hep-th]]. 
[16] J. -W. Chen, J. Deng, H. Dong and Q. Wang, Phys. Rev. D 83, 034031 (2011) [Erratum-ibid. 

D 84, 039902 (2011)] [arXiv:1011.4123 [hcp-ph]]. 
[17] J. Fuini, III, N. S. Dcmir, D. K. Srivastava and S. A. Bass, J. Phys. G G 38, 015004 (2011) 

[arXiv:1008.2306 [nucl-th]]. 
[18] Z. Xu and C. Greiner, Phys. Rev. C 76, 024911 (2007) [hep-ph/0703233]. 
[19] A. El, A. Muronga, Z. Xu and C. Greiner, Nucl. Phys. A 848, 428 (2010) [arXiv: 1007.0705 

[nucl-th]]. 

[20] H. Grad, Commun. Pure Appl. Math. 2 (1949) 331 . 

[21] S. R. de Groot, W. A. van Leeuwen, Ch. G. van Weert, Relativistic Kinetic Theory: 

Principles and Applications, North Holland, Amsterdam, 1980. 
[22] A. Muronga, Phys. Rev. C 76, 014910 (2007) [nucl-th/0611091]. 

[23] E. Byckling, K. Kajantic, Particle Kinematics, Wiley, 1973. 

[24] G. S. Denicol, H. Niemi, E. Molnar and D. H. Rischke, arXiv: 1202.4551 [nucl-th]. 



12 



